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Abstract 

We measure the propagator length in imaginary time quantum mechanics by Monte Carlo 
simulation on a lattice and extract the Hausdorff dimension dn- We find that all local 
potentials fall into the same universality class giving = 2 like the free motion. A 
velocity dependent action (5* oc J dt \ v 1°") in the path integral (e.g. electrons moving in 
solids, or Brueckner's theory of nuclear matter) yields dn = if a > 2 and dn = '2 H 
a < 2. We discuss the relevance of fractal pathes in solid state physics and in QFT, in 
particular for the Wilson loop in QCD. 

PACS index: 03.65.-w, 05.30.-d 
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When considering the geometry of particle propagation in quantum mechanics the suit- 
able language is that of path integrals. Feynman and Hibbs have shown that quantum 
mechanical pathes can be viewed as zig-zag lines, which are no-where differentiable, and 
which exhibit self-similarity when viewed at different length scales. Abbot and Wise 
have shown that quantum mechanical pathes of free motion are fractal curves of Haus- 
dorff dimension two. In imaginary time (Euclidean) quantum mechanics the free motion 
resembles very much the Brownian motion of a classical particle as has been observed by 
Nelson Q]. It turns out that the average path of a classical particle carrying out a Brownian 
motion (steps Ax and At, such that the diffusion coefficient d = ^Ax'^/At is finite) is also 
a fractal curve of Hausdorff dimension two f^, The Euclidean path integral of imaginary 
time quantum mechanics corresponds to a Wiener measure 0. The case of the Euclidean 
free scalar field is discussed by Itzykson and Drouffe [^]. Quite generally, the measure of 
Euclidean path integrals gives the dominant contributions coming from field configurations 
corresponding to pathes of no-where differentiable curves 0. By measuring the length L 
of a curve in terms of an elementary length Ax, the geometrical property of being a fractal 
is captured in the Hausdorff dimension dn, defined by 

L ~A.^o Lo (Ax)'-''" . (1) 

Strictly speaking, this gives a fractal dimension df which in most cases coincides with the 
Hausdorff dimension dn- In quantum mechanics. Ax —>■ requires At — 0, one is at the 
continuum limit. The Hausdorff dimension is a critical exponent. 

In this letter we want to investigate the critical exponent dn of quantum mechanics 
in the presence of local potentials and velocity-dependent potentials in the action of the 
path integral via Monte Carlo simulations on a temporal lattice. In order to be able to do 
a Monte Carlo simulation we consider imaginary (Euclidean) quantum mechanics with a 
negative action 

where we have discretized the time T = tf — ti, At = T/N = e. We define the quantum 
mechanical expectation value of the length for propagating the particle from {xi,ti) to 
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{xf,tf) in D space dimensions via 



^ ^ Jd^xi ■ ■ ■ d^XN-i J2k=o I ^fc+i - I exp[5'[f, e]] 

fe=0 



/ d^xi ■ ■ ■ d^XN-i exp[S'[x, e]] 



It is useful to introduce another quantity related to the length squared 



N-l 



k=0 



XO=Xi,XlSf=Xf 

(3) 
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This can be computed analytically 



< SQf,i{e) >-- 
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In the continuum limit {At — > 0, | Ax |— 0), we expect critical behavior 



(At)-^ ~ (Ax) 



(At) 



-B 



1-dH 



{Axf. 



(6) 
(7) 



In order to determine the Hausdorff dimension, one needs to know the elementary length 
Ax. Prom the action one can estimate this in D = 1 dimension from (Ax/At) ^ ^ 1, 



to give in D dimensions <| Ax |>est~ y2DhAt/m. From that we can estimate < Lf^i > 
and < SQf^i > 



<Lf^i{t)>est = N <\ Ax\>=T^J2Dh/mAt^2DhT/m{<\ Ax\>) ^ , (8) 
<SQf,i{e)>,st = N <{Axf>=2DhT/m. (9) 

Comparing eq.(8) with eq.(6) yields A — \/2 and dn = 2 confirming the known result for 
free motion. Eq.(9) is in agreement (apart from an overall factor 2) with the exact result 
of eq.(5) (for e — > 0) and yields B = C = 0. 

Results of < Lf^i > and < SQf^i > from Monte Carlo simulation are shown in Fig.[l]. 
We have considered in D = 3 space dimensions the free motion as well as motion in the 

presence of a local potential V{x), i.e. the cases: (a) V = (free motion), (b) V = Vqx^ 
(harmonic oscillator), {c) V = Vq/ \ x \ (Coulomb potential). We have plotted log < L > 
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versus log and < SQ > versus A^, where T = NAt is kept fixed. The expectation value 
<| Ax \> is obtained by < L > /N. In Fig.[la,b] we show < L > and < SQ > for the 
free case and find agreement with the estimated scaling behavior given by eq.(6-9). In 
Fig.[lc] we display < L > for the harmonic oscillator. Within numerical errors we find for 
the free motion as well as for both local potentials the same slope. We have extracted the 
Hausdorff dimension dn which in all cases is found to be consistent within numerical errors 
with (i// = 2. Thus in the language of critical phenomena these local potentials fall all 
into the same universality class. For the free motion < SQ > is consistent with the exact 
result from eq.(5). Again the presence of the local potentials does not change the critical 
exponent C which is compatible with zero. 

The above estimate of the Hausdorff dimension hinges upon the presence of the velocity 
squared term in the action. Thus one expects a different fractal dimension if a term with a 
different velocity dependence is present in the action of the path integral. Such a velocity 
dependence occurs in solid state physics for electrons moving in solids ||^. Another example 
in nuclear physics is Brueckner's theory of nuclear matter ||^, using a potential of the form 
V{k) = VQ+V2k'^ + Vik^ . One should note, however, that a non-Gaussian momentum depen- 
dence in the quantum Hamiltonian does not correspond to the same functional form in the 
action of the path integral. Assuming a polynomial momentum dependence, e.g. p^, in the 
Hamiltonian leads in the action to a term which falls off faster than a Gaussian as a func- 
tion of As/ At^/^. This can be parametrized by exp[— (Ax/At^/'^)^^ H — ■ + (Ax/At^/^)^"]. 
If A = 4/3, this leads to a velocity dependence / dtv^^^ in the action. All terms in the 
above parametrization lead to a Ax ~ At behavior different from that of the free motion 
Ax ~ \^At. In order to keep matters simple, we have chosen to consider the following 
parametrization of the action of the path integral 

S = -Jdt ^mv^ + K) K 1° +U{x). (10) 

Here U denotes a local potential. We have seen above that a local potential does not 
change the fractal dimension, so we drop the U term in what follows. Let us estimate the 
Hausdorff dimension for the action S = — J dtV{v). In analogy to the free case we have 
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AtVo{<\ Ax \ > /At)"' oc 1, and thus 

<| Ax \> oc {At)"^ . (11) 

This imphes 

< %,(e) >est cx (At)-^ cx <| Ax |>e~P, (12) 
<SQf,i{e)>est oc (At)"^ oc<|Af|>i^, (13) 

Of course, for a = 2 the scahng behavior agrees with that of the free motion, eqs.(8,9). Prom 
eq.(12), we imply for the Hausdorff dimension dn = ^^ty for the critical exponent C ~ 
• This opens the mathematical possibility of having a very large Hausdorff dimension, if 
a is above but close to unity. However, physically this is not realized, because the presence 
of the kinetic term in the action limits o?// < 2, as will be shown below. 

The results of the Monte Carlo simulations with action given by cq.(lO) corresponding 
to D = 3 dimensions are shown in Figs. [2,3] In Fig. [2] we display < L > and < SQ > 
versus logA^ for a = 1/2. We have extracted the critical exponents du and C as a function 
of a and show a plot in Fig. [3] . We find within numerical errors the following behavior: For 
a <2, dH agrees with the value of the free motion, i.e., dn — 2, while for a > 2, dn agrees 
with the value i.e., the value expected from the velocity dependent potential without 
kinetic energy. Thus acrit = 2 is a critical value of the velocity dependent potential. The 
same kind of behavior is seen for the critical exponent C. 

To conclude let us discuss some perspectives for solid state physics and for QCD. Con- 
sider in solid state physics a semiconductor where conduction properties (due to doping) 
change with position, e.g., like in a transitor. The Hausdorff dimension is sensitive to higher 
order terms in the energy-momentum dispersion relation. Although the determination of 
dn needs a non-local observable, i.e., < L >, the absolute length is not relevant for dn- So 
in principle one could measure < L > on a length scale large to atomic distance, but small 
to the length of variation of macroscopic properties. In such a way, one could determine dH 
locally and study its change as a function of position throughout the semiconductor (this 
would correspond to a velocity dependent potential, where o; is a function of position). 
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In this sense dn can be viewed as one parameter probing the local geometry of quantum 
mechanics. In QCD, a fractal dimension might eventually play a role in a non-local observ- 
able, like the Wilson-loop. The Wilson loop distinguishes between a confined (area law) 
and a non-confined (perimeter law) phase. It also has a quark-quark potential interpreta- 
tion which yields the string tension, being crucial for the baryon spectrum. The potential 
interpretation makes the semiclassical assumption that the quarks propagate along classical 
smooth curves. In quantum mechanics however, one would expect the quarks to propagate 
along zig-zag lines. The path of quark propagation enters in the definition of the Wilson 
loop. Going over from smooth curves, dn — to zig-zag lines, e.g., with dn — which are 
similar to a surface, would most hkely change the behavior of this observable. 

Acknowledgement H.K. is grateful to J. Polonyi and V. Branchina for helpful discus- 
sions. H.K. was supported by NSERC Canada and FCAR Quebec. 
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Figure Caption 

Fig.l (a) Length < L >, and (b) sum of length squares < SQ > for free motion and (c) 
length < L > for harmonic oscillator versus number of time steps N. 

Fig.2 Same as Fig.[l] for velocity dependent action [a — 1/2). 

Fig.3 Critical exponents and C versus potential parameter a for velocity dependent 
action. The full line gives the estimated exponents from eq. (11-13). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9501018vl 
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